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By considering the continuous model for graphene, we analytically study a special gauge field 
for the edge state. The gauge field explains the properties of the edge state such as the existence 
only on the zigzag edge, the partial appearance in the fc-space, and the energy position around 
the Fermi energy. It is demonstrated utilizing the gauge field that the edge state is robust for 
surface reconstruction, and the next nearest-neighbor interaction which breaks the particle-hole 
symmetry stabilizes the edge state. 
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1. Introduction 

A single layer of graphite or graphene is an element 
of various carbon-based materials such as carbon nan- 
otubes.^ These materials are eminently suitable for fu- 
ture nanotechnology and thus fundamental understand- 
ing of the electronic properties is indispensable. The elec- 
tronic property of graphene is mainly determined by the 
TT-electrons around the Fermi energy at the hexagonal 
corners (the K and K' points) of the two-dimensional 
Brillouin zone. It is known that the eigenstates are clas- 
sified into extended states and localized states. The edge 
state existing near the zigzag edge is an example of the 
localized state. ^ The tight-binding (TB) lattice model 
has shown that the edge state exists near the zigzag 
edge while no localized state exists near the armchair 
edge.^'^ Recently, by scanning tunneling microscopy and 
spectroscopy (STM/STS) experiments, Niimi et al.^ and 
Kobayashi et al.,^ respectively, observed the edge state 
near the zigzag edge and the energy position. 

The edge shape determines the boundary condition for 
the wave function. It has been shown in the studies of 
single-walled carbon nanotubes (SWNTs) that a bound- 
ary condition can be changed by an external magnetic 
field and uniform lattice deformations. Using the contin- 
uous model of graphene,^ Ajiki and Ando*" showed that 
a uniform magnetic field parallel to a SWNT axis affects 
the energy band structure of extended states. The mag- 
netic field is expressed by a gauge field which changes 
the boundary condition around a tube axis through the 
Aharonov-Bohm (AB) effect. Zaric et al.^ observed op- 
tical signatures of the AB effect. Kane and Mele^ intro- 
duced another homogeneous gauge field in the contin- 
uous model. The gauge field represents uniform lattice 
deformations such as uniform bend, twist and curvature. 
The gauge field changes the boundary condition around 
a tube axis through a generalized AB effect and makes 



* E-mail: sasaken@imr.tohoku.ac.jp 

t E-mail: murakami@appi.t.u-tokyo. ac.jp 

t E-mail: rsaitoOflex. phys.tolioku.ac.jp 



a mini gap in chiral metallic SWNTs. The curvature- 
induced mini gap was observed by STS experiment by 
Ouyang et al.^^ The perturbations above modify the 
boundary condition for extended states and arc simi- 
lar in that they are represented by gauge fields in the 
continuous model. Is it possible to change the boundary 
condition for the edge state? To answer this question, a 
generalization of the gauge fields to local fields is nec- 
essary since the edge states appear locally around the 
edge. In the previous paper, -'^^ we generalize the gauge 
field introduced by Kane and Mele, to include a local lat- 
tice deformation of graphene. The deformation-induced 
gauge field accounts for the local modulation of the en- 
ergy band gap,^^ which was observed in some pea-pod- 
like structures by Lee et al}^ In this paper, using the 
continuous model, we show that the edge state can be 
formed by a deformation-induced "magnetic" field and 
we clarify why the edge state is formed in the zigzag 
edge, but not in the armchair edge. We explain various 
properties of the edge states, and examine how the edge 
state are affected by perturbations. 

Here we briefly introduce the continuous model. In an 
undeformed graphene, the Hamiltonian for 7r-electrons 
around the Fermi point are given by two-component 
Weyl equation.^ Two components of the wave function 
represent the two sublattices in the unit cell. The Weyl 
equation gives two-dimensional linear k energy dispersion 
relation around the K and K' points in the Brillouin zone. 
The Weyl equation allows us to treat a variety of pertur- 
bations in a unified way for studying transport proper- 
ties, scattering process around impurities, lattice defects 
and topological defects. "'^'^ 

This paper is organized as follows. In § 2, we define a 
gauge field for a locally deformed graphene system. Then 
we solve the continuous model (Weyl equation with the 
gauge field) to obtain the edge state. In § 3, we present 
several applications of the continuous model. The ef- 
fects of a lattice deformation around the edge, the next 
nearest-neighbor interaction, and an external magnetic 
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field on the edge state are examined. In § 4, wc dis- 
cuss and summarize the results. We derive the contin- 
uous model in Appendix A. In Appendix B, we analyze 
a locally deformed graphene system using the TB model. 

2. Continuous Model 

We consider an undcformcd graphene as shown in 
Fig. 1(a). The unit cell consists of A sublattice (•) and 
B sublattice (o), and each sublattice is spanned by prim- 
itive vectors ai and a2. |2a2 — ai| =2£ is the unit length 
along the y-axis. A 7r-electron hops from one site to 
the nearest-neighbor sites with the hopping integral —70 
(« —3 cV) in each hopping process. Let vf{= 7o^/fi) 
be the Fermi velocity and p = {PxiPy) the momentum 
operator. The Hamiltonian for 7r-electrons around the K 
point is then written as^ 





-Px + ipy 



-Px - tPy 




, (1) 



where tr' = {—cTx,<7y) and ai {i = x,y,z) are the 
Pauli matrices. The Hamiltonian operates on a two- 
component wave fimction 'ijj^{r) = '(-(/'^(r), -^g (r))- This 
two-component character is referred to as the pseudo- 
spin; hence, V'a('') ^^'^ ''Psi^) ^'^'^ ^he pseudo-spin up 
and down states, respectively. The energy eigenequation 
HO-V'l(r) = Etp^{r) is the Weyl equation. 



unit cell 



(b) B,{r)^0 




Fig. 1. (a) Lattice structure of graphene. The unit cell con- 
sists of A (closed circle) and B (open circle) sublattices, which 
leads to the pseudo-spin structure of the wave function. Lo- 
cal modulations of the hopping integral is defined by (r) 
(a = 1,2,3). Two examples of local deformation (borders) and 
resultant deformation-induced gauge fields axe shown in (b) and 
(c). In (b), the gauge field (arrows) gives a finite deformation- 
induced magnetic field (flux) illustrated by and (8> while no 
deformation-induced magnetic field is present in (c). 



When the lattice is deformed locally, the nearest- 
neighbor hopping integral depends on the position as 
—70 + ^7o(r) (a = 1, 2, 3) where a denotes three nearest- 
neighbor B sites from an A atom (see Fig. 1(a)). The 
Hamiltonian can be expressed by 



WK=VF<T'-(p-A(r)), 



(2) 



where the vector field A(r) — {Ax{r), Ay{r)) is the 
deformation-induced gauge field. For a weak lattice de- 



formation satisfying |(57a(r)| <C 70, the gauge field is ex- 
pressed by a linear combination of (57a (r) as^^ (see Ap- 
pendix A) 



vpAxir) = (57i(r) - - {5^2{r) + (573(r)) , 



(3) 



Here, we give two examples of local lattice deformation 
on a line (Figs. 1(b) and 1(c)). From now we call this type 
of lattice deformation border. The example of Fig. 1(b) 
is a modification of C-C bond lengths in a row at y = 0, 
namely, (57i(r) ^ at ?/ = and 6^2(1^) = ^73(1") = 0. 
From eq.(3), we obtain Wi?A(j/) = {Sji{y),0). The ex- 
ample of Fig. 1(c) is that 6'ji{r) — and (572(r) = 
^73(1") 7^ at a; = 0. We obtain vfA{x) = (—(572(0;), 0) 
for this case. Figures 1(b) and 1(c) correspond to the 
zigzag edge and the armchair edge if we take Sji = 70 
and 6^2 = Sj^ = 70, respectively. The difference be- 
tween the two gauge fields becomes clear by defining the 
deformation-induced "magnetic" field. 



Mr) 



dAyir) dA^iv) 



dx 



dy 



(4) 



which is perpendicular to the graphene plane B(r) = 
(0, 0, Bz{r)). It should be noted that it is not a real mag- 
netic field. It is called a magnetic field only because it 
enters the Hamiltonian in the similar way as the real 
magnetic field. We have a finite deformation-induced 
magnetic field for Fig. 1(b), and zero magnetic field for 
Fig. 1(c). The gauge field which gives zero magnetic field 
can be removed from eq.(2) by a gauge transformation.^^ 
As shown in Fig. 1(b), the magnetic field appears locally 
around y = 0; Bz{y) is negative at jy < as illustrated 
by ® and is positive at y > as ©. The deformation- 
induced magnetic field will account for the presence of 
edge states not at the armchair edge but at the zigzag 
edge. The importance of the magnetic field can be seen 
by the fact that the energy eigenstate of eq.(2) also sat- 
isfies W^Vl(r) = ^¥|(r): 

4 {(P - A(r))2 + ;iS,(r)a4 Vl(r) = i^VIW- (5) 

has only diagonal-components in the 2x2 matrix. 
One can see that Bzijc) in eq.(5) works as a potential 
and affects the electronic properties. This is because the 
diagonal element is similar to the non-relativistic Hamil- 
tonian Hnr = (l/2m)(p — A(r))^ ± V{r) where the sign 
in front of V{t) corresponds to the pseudo-spin and m is 
a mass. 

We derive the edge states from the Weyl equation with 
the gauge field of Fig. 1(b), corresponding to the zigzag 
edge. Before going into detail, we outline the results here. 
We will show that there are localized states in the energy 
spectrum and the energy dispersion appears as the two 
solid lines at p^, > shown in Fig. 2(b). The velocity 
becomes small with increasing the gauge field and it be- 
comes zero when the gauge field is sufficiently strong. The 
result of the continuous model reproduces the flat band 
of the edge states shown in Fig. 2(e) which is calculated 
from the TB modeP with the zigzag edge. 
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(a) 



E 




rewrite the energy eigenequation as 



Localized state 



(C) C: 



Edge state 
(d) c=l/2 (e) c=l 




Edge state 



Fig. 2. (a) Band structure {px, E) around the K point of the Weyl 
equation without the deformation-induced gauge field and (b) 
that with the gauge field at the border. The shaded regions in 
(a) and (b) represent the spectrum of extended states. The (two- 
dimensional) dispersion relation around the K point for graphene 
without the border gives two cones (representing the linear k 
dispersion) whose apex is the K point as is shown in the inset 
of (a). The band structure of (a) is a one-dimensional projection 
(onto py plane) of the linear k dispersion. In (b), the energy 
dispersion for the localized state is represented by solid lines. 
When the bonds at y = become disconnected and zigzag edges 
appear, the energy eigenvalues of the localized state converge to 
zero (E — > 0) for any Pxi> 0) and forms a flat energy band of 
edge state. Using the TB model, we plot the band structure for 
an undeformed graphene in (c), for a zigzag ribbon in (e), and for 
a graphene with the border (with weakened hopping for the C-C 
bonds at 2/ = 0) in (d). q in (c), (d) and (e) is the momentum 
parallel to the border and g = 27r/3 is the K point. 



We assume that A(r) of Fig. 1(b) is quite locaUzed 
within \y\ < £,g where is a length of the order of lattice 
spacing, namely, A(r) = {Aj;{y),Q) and Aj:{y) = for 
\y\ > £,g in eq.(2). We parameterize the localized energy 
eigenstate as 

^g(r)=^'exp(z^)e-«(^)(^!^J^;), (6) 

where N' is a normalization constant. The pseudo-spin 
modulation is represented by g{y), and the wave vector 
Px is a good quantum number because of the transla- 
tional symmetry along the border. Putting eq.(6) to the 
energy eigenequation, TYkV'bI'") = -^V'B(r), we obtain 



P. - Ax{y) - (G(y) + g{y)) = -—e+'a^y\ 
ay vp 

P. - Axiy) + (G(y) - giy))) = -^e-^^(^) 
ay vp 



(7) 



By summing and subtracting the both sides of eq.(7), we 



Px 



My) 

dGjy) 
dy 



- dg{y) 
dy 



E 

VF 



cosh{2g{y)), 



E 

— sinh(2.g(y)). 

VF 



(8) 



Since we are considering a localized solution around 
the border, we assume G{y) ~ \y\/£, where ^ (> 0) is the 
localization length. When G{y) ^ \y\/£,, the solution of 
the second equation of eq.(8) is given by 



9{y) 



■ i sinh 
i sinh" 



(9) 



The functions G{y) and g(y) are schematically shown 
in Fig. 3(a) and 3(b), respectively. The sign of g{y) 
changes across the border; this sign change means that 
the pseudo-spin flips at the border. The flip is induced 
by the gauge field Ax{y)- To see this, we integrate the 
first equation of eq.(8) from y — — to ^g, and acquire 



'^dy^- 
dy h 



is, 



Ax{y)dv- 



(10) 



We have neglected other terms, since they are propor- 
tional to ^g and become zero in the limit of S^g = 0. By 
putting eq.(9) to eq.(lO), we find 

Ax{y)dy. (11) 



When the right-hand side of eq.(ll) is large, one obtain 
from eq.(9) that g{y) » for y < — and g{y) ^ for 
y>ig- In this case, the localized state is approximately 
a pseudo-spin-up state V'l(r) oc *(1,0) for y < —^g and 
a pseudo-spin-down state V'l(r) oc *(0, 1) for y > £,g. 
Hence, a strong gauge field at the border makes pseudo- 
spin-polarized localized states. Since the polarization of 
the pseudo-spin means that the wave function has am- 
phtude only A (or B) sublattice, this result agrees with 
the result by the TB model for the edge state. ^ 




Fig. 3. (a) The amplitude of the wave function, exp(— G(j/)), of 
a localized state whose localization length is ^. (b) The pseudo- 
spin modulation part, g{y). From eq.(9), g{y) is a constant for 
\y\ ^ Csi S'li'i abruptly changes across the border (y = 0). 



Having described the wave function of the localized 
state, let us now calculate E and ^. To this end, we use 
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the first equation of eq.(8) for |y| > and obtain 

E -px 

Vf 



cosh I 



^ j Ax{y)dy 



(12) 



Moreover, using eq.(ll), we find 

-=Px tanh ( - y Ax{y)dy 



(13) 



In addition to this localized state, there is another lo- 
calized state for the same Px with the same ^ but with 
an opposite sign of E. This results from a particle-hole 
symmetry of the Hamiltonian; azTi-KCTz = —'Hk- By the 
particle-hole symmetry operation, the wave function is 
transformed as tp^^ir) = aztp^ir). 

The normalization condition of the wave function re- 
quires that (, should be positive, which restricts the value 
of Px in eq.(13). Indeed, when Ax{y) is positive, eq.(13) 
means that the localized states appear only at > 
around the K point. This is the reason why the localized 
states appear in the energy spectrum only in one side 
around the K point in Fig. 2(b). On the other hand, the 
Hamiltonian around the K' point is expressed by 



-Hk' =VF<T-{p + A{r)), 



(14) 



where cr = {ax ,o-y). The different signs in front of A(r) in 
eqs.(2) and (14) guarantee the time-reversal symmetry. 
Because of the different signs, a similar argument as we 
used for the K point concludes that the localized state 
appears Px < around the K' point. Furthermore, when 
(1/ft) /f^ Ax{y)dy > 0, £ in eq.(12) becomes zero. The 
zero energy eigenvalue between the K and K' points in 
the band structure corresponds to the flat energy band 
of the edge state. ^ 

When Ax{y) is negative {Ax{y) <C 0), a flat energy 
band appears in the opposite side: < around the K 
point and Px > aroimd the K' point. This condition, 
Ax{y) ^ 0, corresponds to the Klein's edges^^ which are 
obtained by removing A or B sites out of the zigzag edges 
having A or B sites. Calculations on the TB model with 
the Klein's edges^^ also agree with these results obtained 
here. 

There are also extended states in addition to the edge 
states. The energy dispersion relation of extended states 
can be obtained as (E/vp)^ = p'^ + Py, by setting 
G{y) ^ —ipyy/h {\y\ > ^g) in eqs.(6) and (8) where 
Py is a real number. The calculated energy bands are 
given by \E\ > vf\Px\, shown as a shaded region in 
Fig. 2(b). It also agrees with the TB calculation shown 
in Fig. 2(e). One can then regard the localized state as 
a state with a complex wavenumber as follows. Since 
1 — tanh^ .T = 1/cosh^.x, we see that the energy dis- 
persion relation of the localized state between E and ^ 
becomes 



EY 2 
Vf J 



\y\ > i 



(15) 



which is the same as the linear dispersion relation 
{E/vfY = Px~^Py if one replaces Py with i{h/^). 



We have shown that three basic properties of the edge 
states, i.e., the pseudo-spin polarization, the dependence 
on the momentum, and the flat energy band, obtained 
previously by the TB model, ^ can be explained in terms 
of the gauge fleld. In order to quantitatively compare the 
present theory with the TB model, we have performed a 
TB calculation for the geometry of Fig. 1(b) with chang- 
ing 5^fi. Here, we introduce an adiabatic parameter c 
by (57i = c7o, that is, c = and c = 1 correspond 
to no deformation and the zigzag edge, respectively. In 
Figs. 2(c), 2(d) and 2(e), we plot the band structure for 
c = 0, c = 1/2, and c = 1. Comparing these figures with 
the results of the continuous model, one can find a good 
correspondence between the TB model and continuous 
model. Moreover, we analytically find (see Appendix B) 

1^1 _ IpxI 

Vf cosh(— ln(l - 

^ ^ " (16) 

: Px tanh(- ln(l - c)) Oilpl/K), 



+ 0{lpl/h), 



h 



for localized states around the K point. Thus, by com- 
paring eq.(16) with eqs.(12) and (13), we conclude that 
the TB model and the continuous model agree with each 
other near the K point {Pxl/fi "C 1), by the following 
relationship 

Ax{y)dy = -\n{l-c). (17) 

The right-hand side diverges when c — > 1, which repro- 
duces the flat energy band {E 0) in eq.(12) and gives 
^/^ = Px^ in eq.(13). 

When a deformed graphene is located in an external 
magnetic field, the Hamiltonians for the K and K' points 
are given by replacing p in eqs.(2) and (14) with p — 
eA'''^(r) as 

Wk = vf<t' ■ (p - eA^-(r) - A(r)), 

WK'=^^FO--(p-eA<^-(r)+A(r)), ^^^^ 

where e denotes the electron charge and A''™(r) is 
the gauge fleld for the external magnetic field. The z- 
component of the external magnetic fleld is given by 



sr(r) 



5A7(r) aA°'»(r) 



dx 



dy 



(19) 



The same signs in front of eA°™(r) in eq.(18) reflect the 
violation of time-reversal symmetry. Due to the similarity 
between the external magnetic field and the deformation- 
induced magnetic field in eq. (18), one may expect that 
the edge states can be induced also by A®™(r). However, 
as we will see in § 3, the deformation-induced magnetic 
fleld for the edge state corresponds approximately to one 
flux quantum in a hexagonal unit cell, and is thus ex- 
tremely strong (^ 10^ T), compared with the external 
magnetic fleld. 

3. Application of Continuous Model 

In this section, we apply the continuous model to the 
following three examples. 

First, we consider an effect of surface reconstruction 
(SR) around the zigzag edge on the edge state. Since the 
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SR is an additional lattice deformation around the edge, 
it can be expressed by an additional gauge field in the 
continuous model. Let 5A^^(r) = {5A^, 5A^^) be the 
gauge field for SR. The perturbation is then represented 
by 

+ i5Al^\ 

Af )■ 

(20) 

To calculate the energy shift by 5'H^^{t), we recall 
that the edge state (V'ij^o) ^ pseudo-spin polarized 
state. Since STi^^ir) only mixes the wave function on 
A and B sublattices to each other, it is shown within 
the first-order perturbation theory that (57i^^(r) does 
not shift the energy of the edge state, namely, /S.E^^ = 
{il)f\5H^^\-4)f) « 0. This resuh is consistent with the 
result of Fujita et al.,^^ who found that additional lat- 
tice distortion has little effect on the energies of the edge 
states. 

Second, we consider next nearest-neighbor (nnn) hop- 
ping. The nnn interaction decreases the energies of (i.e., 
stabilizes) the edge states, as shown by the TB calcula- 
tion.^® It accounts for the STM/STS measurements,"''^ 
in which the edge state appears about 20 mcV below the 
Fermi energy. We note that the importance of the nnn on 
a localized state around a defect is pointed out by Peres 
et alP 

The reason for the stabilization is that the nnn inter- 
action mixes the wave function on the same sublattices, 
namely, it gives diagonal terms in eq.(2). Here, we pro- 
pose a simple way to include the nnn hopping process 
into the continuous model, and then calculate the energy 
shift of the edge states. Let — 7„ (« — O.I70) denote the 
nnn hopping integral. For an undeformed graphene, the 
nnn interaction can be included by adding — 7n(V^)^P^ 
to eq.(l) as^*^ 



where (o'z(r)) = {^p^ {v))'^ a z'4>^ {r) is the pseudo-spin 
density. It follows from eqs.(6), (9), (11) and (17) that 



vpcr • P - 7n I ^ 



P • 



(21) 



The border is incorporated by replacing p with p — A(r) 
in eq.(21). As a result, we can write the total Hamiltonian 
as 

HK=HK-7n(0 (P-A(r))2. (22) 

It is noted that this replacement is not always correct 
in general because the modulation of the nnn hopping 
integral is generally independent of the nearest neighbor 
one. In the present case, however, we adopt this replace- 
ment as a simplest approximation. Using (Hk/vf)'^ = 
{p — A{r))'^ + hBzir)az (see eq.(5)) ineq.(22), we obtain 

^)'+7n^S,(r)a,. (23) 

Having calculated the energy spectrum of TCk, we can 
evaluate the energy shift due to the nnn interaction using 
eq.(23). Because the edge state satisfies Hk^^{tc) = 
{E = 0), the first-order energy shift is given by 



Hk = Hk — In 



I 



1 +^ 

1 — ±M. 



{y < 0), 
(y > 0), 



(25) 



/ / Bzir){az{r))d\ (24) 



which is nonzero only around the border, reflecting the 
pseudo-spin polarization. Since Bz{y) is confined near 
the border as shown in Fig. 1(b), we assume Bz{y) = 
{h/£)S{y - 0+) - {h/i)6{y - 0_). Putting eq.(25) into 
eq.(24), we obtain 

Ai;ji"" = -27„^, (26) 

which reproduces the previous result of the TB model 
for the edge state. -'^ 

An important feature of the nnn interaction is that 
the energy shift is always negative and therefore it always 
stabilizes the edge state. The stability comes from the in- 
teraction between polarization of the pseudo-spin and the 
deformation-induced magnetic field. Another important 
point is that the nnn interaction breaks the particle-hole 
symmetry. As discussed previously, the edge state ap- 
pears in pair (ijj^ and (TziJj^) at E = and their pseudo- 
spin densities are the same. Thus, AE^'^ is the same for 
the two states, and they remain degenerate at the energy 
^^nnn Q) though the particlc-holc symmetry is lost. 
It is clear in eq.(23) that the particle- hole symmetry for 
the edge state is broken locally at the edge by the last 
term of the right-hand side. 

Third, we consider an external magnetic field. Since 
the sign in front of A'^™(r) is the same for the K and K' 
points as shown in eq.(18), the energy shift by a magnetic 
field becomes 

AE«- = ^r.jjj eBr{r){az{r))<fr. (27) 

Thus, i?|™(r) can shift the energy spectrum of the edge 
state. When B|™(r) is uniform, we estimate the energy 
shift as AE^"^ w 7n(^/^B)^, where £b is the magnetic 
length defined by = \/fi/eB^™. Here, we assume that 
eq.(25) is still valid even in the presence of a weak uni- 
form magnetic field. The numerical value of Ib {£) is 
about 25/^5°™ nm (0.2 nm) where the magnetic field 
is measured in the unit of Tesla. For example, Bf^ « 
10 Tesla gives a tiny shift A£'«™ w 0.2 meV. This value 
is one-order smaller than the Zceman splitting. Thus, 
compared with the deformation-induced magnetic field, 
external magnetic field has little effect on the edge state. 

4. Discussion and Summary 

By considering the edge state using the continuous 
model, we found that the deformation-induced gauge 
field (A(r)) and magnetic field (^^(r)) explain basic 
properties of the edge state. It is summarized as follows: 

(1) The gauge field can generate the edge state in energy 
spectrum, depending on the gauge field direction. Let e|| 
the unit vector along the border and A||(r) = A(r) • ey, 
localized states (the edge state) appears if the gauge field 
has a component parallel to the border: Ay (r) ^ 0. The 
edge states are pseudo-spin polarized. 

(2) The direction of the gauge field, namely, An (r) > 
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or j4||(r) < 0, is vital for the edge states, as it determines 
the energy dispersion and wave vectors which allow the 
edge states. 

(3) The flat energy band of edge states at zero energy re- 
sults from a divergence of the gauge field; ^||(r) ±oo. 

(4) The nnn interaction lowers the energy of the edge 
state below the Fermi energy. This is because the nnn 
induces a linear coupling between the pseudo-spin and 

(5) Bz (r) for the edge state corresponds to an enormous 
magnetic field ~ lO^T at the border. Thus a uniform ex- 
ternal magnetic field has little effect on the edge state, 
compared with the deformation. 

Here, we discuss impurity effects on the edge states, 
as impurity scattering may give rise to the Anderson lo- 
calization for low-dimensional systems. However, it has 
been proposed that the Anderson localization does not 
occur in graphene, since graphene has a remarkable prop- 
erty, leading to the absence of back-scattering.^^ Because 
a rotation of the wave vector around the Fermi point 
gives an additional phase shift of tt (Berry's phase), a 
time-reversal pair of scattered waves cancel with each 
other. However, in the presence of deformations such 
as a vacancy, the wave function gets an extra phase shift 
through the AB effect of the deformation-induced gauge 
field. The extra phase makes an asymmetry between 
time-reversal pair of scattering waves and thus the back- 
scattering which contributes to the localization will be 
recovered. 

Finally, we point out a possible experiment to confirm 
the results of the present theory. A direct way of mak- 
ing localized states is by applying a stress locally to an 
undeformed graphene, thereby inducing a deformation- 
induced magnetic field. The present theory predicts that 
a formation of localized states depends sensitively on the 
direction of the stress, and such localized states can be 
probed by a peak in the local density of states in STS 
measurement. 

In summary, we have formulated the edge state as a 
localized sohition of the continuous model for a deformed 
graphene. We found that the deformation-induced gauge 
field at the border is relevant to the presence of the edge 
state. The gauge field can reproduce basic properties of 
the edge state; the pseudo-spin polarization, the disper- 
sion, and the flat energy band. The gauge fleld is induced 
by local lattice deformation in general and also can be 
simulated by an external magnetic field. Compared with 
an external magnetic field, local lattice deformations can 
generate a strong field of order of 10^ T. The gauge field 
description was extended to include the SR and the nnn 
interaction. In terms of the gauge field, we showed that 
the SR does not shift the energy of the edge state and 
the nnn stabilizes the edge state energy. 
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Appendix A: Gauge Field for Weak Lattice De- 
formation 

First, we derive the gauge field quoted in eq.(3). It can 
be obtained by considering the matrix element of the de- 
formed Hamiltonian between the Bloch wave functions. 
The deformed Hamiltonian is defined by 

Wdeforxn^^ ^ d^TiKcf+J ct + h.C, (A-1) 

iGA a=l,2,3 

where Cj and c| are the canonical annihilation-creation 
operators of the electrons at site i. The off-diagonal ma- 
trix element of Hdeform is given by 

(f^+^''|Wdeform|*|) 

= 7^E E ^7a(rO/a(k)e-^^''-% (A.2) 

iGAa=l,2,3 

where /a(k) = e*'^'^" and Rq (a ~ 1,2,3) are vectors 
pointing to the nearest-neighbor B sites from an A site. 
Here, \^^) is the Bloch wave function defined as 

|4'^) = ^^e^''-'cl|0), {s = A,B) (A.3) 

where denotes the number of graphene unit cells. 

By expanding /a(k) in eq.(A-2) around hp (wave 
vector of the K point), we obtain /a(k) = fai^p) + 
i/o(ki;')(k — ki?) • Rq -I- • • • . The second term and the 
higher order corrections can be ignored if we consider 
the low-energy properties near the Fermi points of a de- 
formed graphene satisfying |(57a(r)| <C 7o- Then the dom- 
inant contribution is given by 

l^E E ^7a(rO/a(kF)e-'^''--. (A.4) 

" ieAa=l,2,3 

The K point satisfies kp ■ ai = —Att/S and kp ■ (a2 — 

ai/2) = (mod 2n), and therefore we obtain fiCkp) = 1, 
./2(kF) = e+*"r and /3(kir) = e~*"?. Substituting these 
into eq.(A-4), we see by the definition of A(r) (eq.(3)) 

^^Y.iMri)+iAy{r^)}e-^''^-K (A-S) 

By the same calculation for (^^'''^''IWdeforml'S'A)) ob- 
tain the matrix for Tidcform as —vpa' ■ A(r). Similarly, by 
expanding the nearest-neighbor Hamiltonian for an un- 
deformed graphene around the K point, we see vpa' ■ p. 
Thus, we obtain Vpa' • (p — A), which is eq.(2). 

Appendix B: Derivation of eq.(16) 

Next, we consider a general border using the parameter 
c to obtain eq.(16). We write the Hamiltonian Ho for an 
undeformed periodic graphene shown in Fig. B-1 as, 

Wo = -70 E E(c?+a)^^i' + l^-^- 

a=l,2,3 ieA 

The Hamiltonian for graphene with zigzag edge is de- 
noted by Tizigzag, and we define Tiborder by TLq - Tizigzag- 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 7 



'^border describes the hopping across the border. Prom 
the definition, Tiborder is given by 

^border = "70 ^^(C/, j=jv)^cf_^=o + ^-C., {B-2) 
I 

where c/^j and c\ j are annihilation-creation operators 
of TT-electron at the edge site, (/, J) ( J = or N). In the 
following, we consider 7i(c) = Hq — cTibordcr where \c\ <^ 
1 corresponds to a weak deformation and c ^ 1 leads 
to Wzigzag, i-e., zigzag edge formation. In addition to the 
zigzag cdge^ for c = 1, Ti{c) describes the Klein's edge^^ 
when c —oo. When c — oo, an electron cannot 
come to the border sites, and the system reduces to a 
graphite ribbon with the Klein's edges which are made 
by removing the border sites from the zigzag edges. 

Before we examine localized state of 7i(c) analytically, 
let us show the numerical result of the energy band struc- 
ture around the K point including both localized state 
and extended state. In Figs. 2(c), 2(d) and 2(c), we plot 
the energy spectrum of H(c) for = 10 with c = 0, 1/2, 
and 1. Because Hbordor does not break the translational 
symmetry along the border, the wave vector, g = k • ai, 
remains to be a good quantum number. We consider the 
spectrum near q = 2ir/3 (the K point). In the energy 
spectrum of H{0) (Fig. 2(c)), the valence and conduction 
bands touch a,t q = 2n/3. Only extended states exist in 
the system and they are doubly degenerate, which are 
supported by the (inversion) symmetry under py ^ —py 
where Py denotes the momentum perpendicular to the 
border. These degeneracies disappear as a result of the 
deformation, and the spacing of energy bands for 7i(l) 
becomes about one half of that for Ti{0). Apart from 
the change of energy spacing, we observe no significant 
changes in the band structure of the extended states. 



J = 



Border 




J = N 



Fig. B-1. Geometry of graphene with the zigzag border. When 
c = 1, the system becomes a graphene with zigzag edges since 
the hopping integral across the border is —(1 — c)7o. 



On the other hand, the localized states have several 
notable characteristics. First, localized state appears for 
27r/3 < < TT as we change c from to 1, i.e., there 
is an asymmetry about the K point in the momentum 
axis. Second, the two localized states gradually depart 
from the valence and conduction bands respectively as 
c is increased (Fig. 2(d)). Third, those localized states 
merge when c = 1 to become the flat bands with zero 
energy (Fig. 2(e)). In the following, we will analyze these 
features by analytically calculating the energy eigenvalue 
and localization length of localized state. 



An energy eigenstate can be represented as 

^i,(/,J) = iV'exp(zg/) ("^^l;^! 



(B-3) 



The energy eigenequation Ti.{c)ijjE = EipE is written as 

eMJ + 1) = -MJ) - QM-J + 1), 



(B-4) 



where Q = 2 cos(g/2) and e (= E/-jo) is the energy eigen- 
value normalized by the hopping integral. The eigenstate 
is expressed as 

S (^^^ . ' ^ 



^a(J) 




{J = 0,...,N) 



(B-5) 



where Sj = sin( J^)/ sin Here is a wavenumber in 
the direction of J and satisfies 



(3^ + 2Qcos(/>+l. 



(B-6) 



The derivation of eq. (B-5) is as follows: first, it is clear 
from eq.(B-4) that eq.(B-5) holds for J = and J = 1. 
Second, let us rewrite eq.(B-4) using eq.(B-6) as 



Q + 2 cos t 



Qj \MJ) 



(B-7) 



Equation (B-5) can then be proved by induction. 

The wave function is given by three parameters c, q 
and (j). The parameter is determined by the boundary 
condition: the energy eigencqTiation at border sites, J = 

(B sites) and J = iV (A sites). 



e,^B(0) - -(1 - c)MN) - QMO), 
ecPAiN) = -{l- c)MO) - QMN). 
Equation (B-8) can be rewritten as 




(B-8) 



• (B-9) 



From eqs.(B-5) and (B-9), we obtain two equations re- 
lating {(f>A{0),MO)) to (0A(iV),<^B(iV)). Then, by elim- 
inating (/>A(iV) and 4>b{N) in these equations, we obtain 
the following constraint equation for 0, 



det 



where 



Sn+1 \ _|_ Q 
Sn ) ~'~ An ''^'^J An ' 



Sn- 



Sn 




{1-c)Sn, A: 



{1-cf 



(B-11) 



Equation (B-10) has 2N solutions of (j) for given N. From 
the numerical calculation, we find that for > 1, 
namely, for 0< q < 27r/3or47r/3 < q < 2tt, the equa- 
tion has 2A^ real value solutions of (p, and for \Q\ < 1 
(27r/3 < q < 47r/3), it has 2(A'' - 1) real value solu- 
tions (extended states) and 2 complex value solutions 
(localized states). As for localized state, ^ can be pa- 
rameterized as (j) = i£/^ or (j) = w + i£/^ depending 
on — 1 < Q < and < Q < 1, respectively. Here, 
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^ is the localization length and thus when ^ ^ oo the 
state becomes an extended state. From eq.(B-6), we have 

= (Qi 1)^ for (/) = and tt. Thus, one can understand 
that there are two points in {q, (j)) plane satisfying e = 0, 
where a localized state turns into an extended state. ^° 
One point is located at (27r/3,7r) {Q = 1) which is the 
K point and the other is located at (47r/3,0) {Q = —1) 
which is the K' point. 

When N ^ 1,'by solving eq.(B-lO) for localized state, 
we acquire 



£ , / c(2 - c) + Vc2(2 - c)2 + 4g2(i _ c)2 \ 

r^H ^\ J ' 

(B-12) 

where • • • denotes correction of 0{e^^^^^) and therefore 
it is negligible when N ^ ^/£. To obtain eq.(B-12), we 
neglected A]^^ in eq.(B-lO) since is large as 0{e^^/^). 
In this approximation, we get Q{Q -\- e~^^^) + A{Q + 
e'^/i) = 2e2 for Q < 0. Using eqs.(B-6) and (B-11), we 
rewrite this condition as 
c(c - 2) 



Q 



{i-cY 



0, 



(B-13) 



which gives eq.(B-12). The energy eigenvalue is obtained 
by inserting eq.(B-12) into eq.(B-6). It is to be noted that 

not only c = but also c = 2 give no localized state since 
eq.(B-13) gives = in these cases. For c = 2, the hop- 
ping integral across the border is +70, while other bonds 
arc —70. By a gauge transformation, this difference in 
sign of the hopping integral on the border can be ab- 
sorbed into an AB phase tt around the nanotube (to the 
y-dircction in Fig. B-1). Thus the system is equivalent 
to a SWNT without border, with a tt flux (a half of the 
flux quantum) parallel to the NT axis. Because it has no 
border, it has no edge state. 

Next, to understand E and ^ for localized state near 
the K point, we define the wave vector fc^, measured from 
K point as q = 271/3 + \ai\kr^. Because Q = 2cos{q/2), 
we have Q=l-ik^ + {Ik^f /% + 0{{ik^f). Then, from 
eq.(B-6) with (p = w + i£/^, we obtain 



(I' 



(B.14) 



where • • • indicates higher order corrections such as 
iPx/fi- The first two dominant terms in the right-hand 
side reproduce eq.(15). By expanding eq.(B-12) in terms 
oi Px{= Kkx), we obtain eq.(16). 

Finally, we check whether E and ^ of the continuous 
model reproduce the results of the TB model for local- 
ized states satisfying k^l ~ 1. First, using eqs.(B-6) and 
(B-12), we plot E and ^ for g = 137r/18, 77r/9 and Stt/Q 
as solid curves in Figs. B-2(a) and B-2(b), respectively. 
In Fig. B-2(a), E decreases slowly at small c. However, 
every curve decreases almost linearly as a function of 
c for large c and converges to zero. In Fig. B-2(b) we 
see that ^ decreases quickly with increasing c and con- 
verge at = -l/ln|Q| (eq.(B-12) with c = 1). The 
convergence agrees with the previously published litera- 
tures. 2' The dashed lines in Figs. B-2(a) and B-2(b) arc 
E and ^ derived from the continuous model. The wave 



vectors q = 137r/18, 77r/9 and Stt/Q, which we used in the 
plot for the TB model, correspond to \aLi\kx = 7r/18, 7r/9 
and 27r/9, in the continuous model. Using eqs.(12) and 
(17), we plot E/"fD (= e) as dashed curves in Fig. B-2(a). 
One can see that the solid and dashed curves are almost 
identical for q = 137r/18. However, their difference is vis- 
ible for q = 77r/9 and not negligibly small for q = 87r/9. 
Using eqs.(13) and (17), we plot ^/i as the dashed curves 
in Fig. B-2(b). The continuous and the TB models agree 
with each other with c near zero. For general c and q, 
while the trend of £, agrees for the two models, their dif- 
ference becomes large as the state is located apart from 
the K point. The discrepancy between the TB and con- 
tinuous models comes from the assumption of the linear 
k energy dispersion relation of the Weyl equation. The 
difference increases as the state is located apart from the 
Fermi point, which represents the deviation of the energy 
band from the linear energy dispersion. The energy dis- 
persion relation may be improved to reproduce the TB 
model by adding some higher order terms to the Weyl 
equation. 
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Fig. B-2. (a) Energy eigenvalue -E/70 and (b) localization length 
as a. function of c. The solid and dashed curves represent the 
TB model and the continuous model, respectively. Three wave 
vectors q = 137r/18, 7n/9 and 87r/9 are chosen as examples. 

When c = 1, each curve in (a) converges to zero, which agrees 
with the flat band of edge states, and each solid curve in (b) 
reproduces previous results of the localization length.^' 
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